Constraints on measurement-based quantum computation in effective cluster states 
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The aim of this work is to study the physical properties of a one-way quantum computer in an 
effective low-energy cluster state. We calculate the optimal working conditions as a function of the 
temperature and of the system parameters. The central result of our work is that any effective 
cluster state implemented in a perturbative framework is fragile against special kinds of external 
perturbations. Qualitative aspects of our work are important for any implementation of effective 
low-energy models containing strong multi-site interactions. 
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A cluster state |^cs) is a quantum state defined on 
some lattice (we focus on a square lattice) of qubits, 
which fulfills the following eigenvalue equations 

K a := X a (g) Z b , K a l^cs) = l^cs) , (1) 
ber(a) 

with X a and Z^ being Pauli operators acting on qubits 
a, 6, and T(a) is the set of nearest neighbors of site a. 
So called one-way quantum computers (1-WQC) perform 
universal quantum computations just by one-qubit mea- 
surements on a cluster state [1]. While this saves the 
need to apply unitary transformations to the quantum 
register, it requires to reliably prepare a cluster state. 
One possibility is to cool the Hamiltonian 

H cl :=-^K a (2) 

aec 

into its non-degenerate ground state, which is by defini- 
tion the cluster state. A direct implementation of H c \ is 
not realistic, as it contains multi-qubit interactions which 
are not realized in nature. 

This suggests to search for a more realistic Hamilto- 
nian with only two-qubit interactions having the same 
non-degenerate ground state. But such an Hamiltonian 
does not exist [2], so one is limited to approximate the 
cluster state. One approach is to use ancillary qubits to 
effectively mediate the many-qubit interactions. Unfor- 
tunately, this is not of practical use, since the necessary 
precision scales with the system size [3]. 

Alternatively, one can implement H c \ as an effective 
low-energy model of a realistic Hamiltonian contain- 
ing only two-qubit interactions [4, 5]. Two questions 
arise naturally: What are the optimal working condi- 
tions to perform measurement-based quantum computa- 
tion (MBQC) in an effective cluster state? How robust 
are effective cluster states with respect to external per- 
turbations? In the following, we show that any such effec- 
tive cluster state, which is implemented in a perturbative 
framework, is strongly affected by external perturbations. 

Model — We replace each qubit on a square lattice 
by 4 physical qubits and we encode the logical clus- 
ter qubit into the subspace defined by the projector 




Figure 1. a) Physical qubits (black dots) on a CaVO lattice. 
The 4 physical qubits of a lattice site fi G C (gray circles) 
are named with the double indices (/i, 1), (/x, 4). The two 
physical qubits of a bond m G M (red lines) are called m(l) 
and m(2). To a physical qubit (/x,i), the neighboring qubit 
on the bond is i). The ZZ-interactions of H are colored 
blue, b) Energy spectrum of H l ° c (Eq. 3). 

P := |0io g ) (0000| + |li g) (Hll| [4]. The Hamiltonian 
H := gHo + \ XZ V is then defined by: 

9Ho := -9 Z M ® Z ^,j) > 

4 

KzV := —X xz ^2 X^O'*) ® ^£0M) • 

/i^C i—l 

The symbol i j means, that the lattice sites i and 
j are related in some connected graph structure. Other 
notations are illustrated in Fig. la. The ground-state 
space of Hq is the space of all logical qubits. 

It is possible to solve this model exactly by trans- 
forming H into the base S loc with the base transforma- 
tion CZm (controlled-Z operators on every bond) [5]. 
Using (CZ M )X M <g> (CZm) = X (/M) <g> 

and [CZm, Ho] = 0, the transformed Hamiltonian reads 
H Xoc = Y.^ C H l ° c with 

4 

H ]? c : = -9 ^2 Z M ® z (ti,j) - ^xz ^2 X M - ( 3 ) 

The H l ° c are local transverse field Ising models (TFIM) 
on each lattice site \i which can be solved by exact diago- 
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nalization. If is the 16 x 16 matrix of the 16 eigenvec- 
tors |0)^ ... 1 15)^ of i^ oc , the diagonal form of H in the 
basis £ dia s is: H dia % := [<g) MGjC i^]-ff loc [® MeC (see 
Fig. lb). Most importantly, the gap between the 
(unique) ground state |^#) and the first excited state 
arises perturbatively in order 4 in X xz /g [5, 6]. 

It is useful to generalize the cluster stabilizers K a into 
the physical space 

4 
i=l 

The if M can be also transformed into the basis H loc : 

K l ™ = (CZ M )K fl (CZ M ) = ®ti X M- In this basis ^ 
is easy to show, that the K l ° c commute with each other 
and with H loc . Consequently, the eigenvalues ±1 for each 
stabilizer are conserved quantities. For \ xz — >> oo, the 
ground state of H loc is the polarized state which is eigen- 
vector to all K l ° c with eigenvalue +1. The ground state 
of H in the limit X xz — >• is therefore the cluster state of 
the logical qubits. 

The two low-energy states {|0) , of H l ° c repre- 
sent an effective qubit on the according lattice site \i (see 
Fig. lb). We can therefore derive an effective low-energy 
model in the space S of all effective qubits 

tf dias l 5 = -^E^ = -^E^ ias l 5 > w 

where being the Pauli Z-operator acting on the effec- 
tive qubit on the lattice site fi and Kf** = R^K^R^ 
It follows, that the effective low-energy approximation of 
H in the limit \ xz — » is the cluster state Hamiltonian 
H c \ of the logical qubits. 

Fidelity — The usability of H for quantum compu- 
tations depends on the question how "well" the logical 
cluster state is approximated by \i/jh)- This can be quan- 
tified by the fidelity F = | (V^l^cs) | 2 of two states and 
its generalization for density operators F = (^csIpI^cs) 
[7] . The fidelity translates to the "success probability" of 
a MBQC using \^h) as a resource state [5]. 

As shown in [8], for a Hamiltonian Hn(X) with 
control parameter A, size TV, and two ground states 
tMA^tMA') it is: lim iV ^ 00 F(^iv(A),^iv(A / )) = d N 
with d G [0, 1] being constant. The "fidelity per site" 
d := limAT^oo a/^XV^A), ip N {\')) is therefore intensive. 
Since d < 1 for X xz > 0, the fidelity of the logical clus- 
ter state with \i/jh) vanishes for N — >• oo. This questions 
the usability for large systems. However, the concept can 
still be applied by quantum error correction techniques 
[5]. In this context the value d translates to the success 
probability per measurement which must be large enough 
to fulfill the threshold theorem [9]. 

We calculate the fidelity F(|^cs) ? p) = (^csIpI^cs) of 
the logical cluster state |^cs) with the canonical density 



operator p := \e~^ H = ^Y,i e ~^ E ^ i} \^%) (^1- Here 
Z = Tr(e~P H ) denotes the partition function, /3 = j^, 
and one finds 

F(\i/>cs) , P) = F{\+) c , pt) = F(\+) fl , p l p N =: d N . 

with := 1/v^dOiog}^ + |lio g ) M )- Consequently, it is 
sufficient to study a single lattice site (we omit index p) 

15 

d = (+|i?t p dia gi? | +) = J- e-^i«> I (i\R\+) | 2 , 

i=0 

with Z = ^2 i e~^ E ^ . Next, we approxi- 
mate e-^i*>|(i| J R|+)| 2 « Wi ^ and 
Z « e -P E \v + e ~ f3E \^ which is justified by the 
following observations: 

(a) For MBQC, we have to choose the temperature low 
enough so that even the first excited state plays a minor 
role. Due to the exponential scaling factor we can omit 
all contributions of high-energy states. 

(b) Due to the orthogonality of the vectors it is 
| (i\R\+) | 2 < 1 - | (0|i?|+) | 2 for all i G {1, ..15}. For not 
too large X xz we expect | (0|i?|+) | 2 < 1, so the contribu- 
tions of the other states are small. 

(c) For the first excited state it is | (1\R\+) | 2 = 0. This 
is proven by the fact, that R |+) and |1) do not have the 
same conserved eigenvalue of the if -operator. 

The resulting fidelity per site d 

^ 1+e U lW+)' 2 ( 5 ) 

is shown in Fig. 2 for different T. For finite T and 
for small \ xz , the fidelity is dominated by thermal fluc- 
tuations. For large X xz the curve follows the zero- 
temperature ground-state fidelity. In between, there is 
a trade off between both effects. If one assumes an 
error correction algorithm for a 1-WQC with a sim- 
ple error model of Pauli errors [10] as it is given in 
[11] with an error threshold of 1.4% (d > 0.986), 
then the maximum T max where this threshold holds is 
T max = 2.18 10~ 4 g/k B . It is reached for \° x f = 0.222 g. 

The Hamiltonian H can therefore be used as a 1-WQC 
under conditions that in principle could be prepared in 
a laboratory. Next, we determine the robustness of such 
an effective cluster state against external perturbations. 

Z -field — First, we consider the presence of an external 
field in Z-direction: 

4 

H z := gH + \ XZ V - h z ^ 

i=l 

Let us formulate H z using the basis £ diag limited to S. 
Since the perturbation commutes with CZm, ° ne nas 

A F 

Hf^\ s = -—Y / ^-^h z c(\ xz / 9 )Y / ^ , 
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— T=0 

_ T=0.0002[g/k B ] 
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T=0.005[g/k! 

— T=0.025[g/k! 
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Figure 2. Fidelity per site d (Eq. 5) in dependence of X xz for 
different temperatures T. Inset: Contribution c(X xz /g) of a 
physical Z-operator to the effective X-operator. 



where c(X xz /g) G R can be read easily from the matrix 
RpZfa^K^. One finds lim^^o c = 1 and the space <S is 
decoupled from the high-energy space for this limit. The 
value c(X xz /g) is plotted in Fig. 2. We stress that the 
scale AE is of order 4 in X xz , while the scale of the X- 
field is proportional to h z . One therefore expects a polar- 
ization of the ground state for very small ratios h z /X xz . 

This is confirmed by solving the Hamiltonian H z ex- 
actly in the basis S loc 



H 



loc 




(am) 



which is still a sum of local terms H l °^. 



The fidelity 

per site of the ground state is calculated as in the un- 
perturbed case using the eigenvectors of H l °^. It can be 
seen in Fig. 3a that already very small ratios h z /X xz have 
a significant impact on d. For the above example, one 
finds the upper bound h™ ax = 1.5210 -5 # satisfying the 
threshold d > 0.986 at T — 0. Thermal fluctuations play 
only a minor role, because the gap is strongly increased 
by the external field [6]. 

ZZ -coupling — Second, we consider the effect of addi- 
tional Ising ZZ-couplings on the bonds m G M 



H zz := gHo + X XZ V — X z 



■■ £ - 



i(l)^m(2)- 



Now we formulate the Hamiltonian using the basis S diag 
limited to S 



meM 



where again the space S is decoupled from the high- 
energy space for the limit X xz — » 0. This Hamilto- 
nian represents a TFIM on the square lattice. For 
this model a quantum phase transition takes place at 

2A I c 2 

— zz ^ t — = 0.3285 separating an ordered from a dis- 
ordered phase [12]. The gap closes at the critical point 




Figure 3. a) Ground-state fidelity per site d — \ (0 z \R z \-\-) | 2 
of H z as a function of X xz and h z for T — 0. b) Fidelity 
per site d (Eq. 8) of H zz in dependence of X xz and X zz for 
T = 0.001g/k B - 



changing the ground state significantly, so the system is 
not useful for MBQC anymore. The energy AE is a 4th- 
order term in X XZ1 while the Ising part is of the order 
X zz . Very small values X zz /X xz are therefore sufficient 
to destroy the cluster phase. The estimated critical line 
^zz\crit is shown in Fig. 3b as a function of X xz . 

To approximately calculate the fidelity per site for fi- 
nite temperatures we use the analog of Eq. 5: 



1 



1 



^ff^e-^dk 



d(\^ Hzz \ s )A^cs\s)), (6) 



We note that point (c) is no longer valid as K is no longer 
conserved, but the use of Eq. 5 is still justified by (a) and 

(b). The dispersion uo(k) w \J (AE ZZ ) 2 + (v • |£|) 2 of the 
first excited mode is taken into account by this equa- 
tion, while correlated excitations are neglected [13]. The 
energy gap AE ZZ of the TFIM is calculated by a dlog- 
Pade [6,6] approximation of its order 13 series expansion 
[12] and v = 0.99 AE/2 [14] is the spin wave velocity at 
the critical point. 
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We now transform l^cs) into the effective basis 



®|Om><o m I^ 



(7) 



(using (11^1+^) = 0) such that Eq. 6 reads 
l(0|i?|+)| 2 



d; 



d 



TFIM 



(8) 



where ^tfim := d{\^ Hzz \ s ) , ®^ eC |0 M )) corresponds to 
the ground-state fidelity per site of a TFIM with the 
polarized state. We have calculated it as a high-order 
series expansion about the high-field limit 



^TFIM 



where A := 



1 



93 
256' 



A 



2961 
2048' 



A b 



243005 



32768 



A 



812949139 
" 18874368 ' 



10 



17716040461601 
65229815808 



12 



2X Z 



AE 



[6]. 



The fidelity per site of H zz \ e s is plotted in Fig. 3b. 
One sees, that close to the critical point the fidelity drops 
due to quantum and thermal fluctuations. For the above 
example, one finds the upper bound A™ ax = 8.33 10 -5 g 
satisfying the threshold d > 0.986 at T = 0. 

We additionally calculated the energy gap and the 
ground state fidelity of the full Hamiltonian H zz as series 
expansions. The high-energy contributions turned out to 
be negligible corrections to the low-energy results [6]. 

Conclusions — We have seen that the effective cluster 
state of H could be used as a 1-WQC under conditions 
that in principle can be prepared in a laboratory. But to 
be of practical use, effective cluster states must be also 
robust against additional perturbations. 

We have shown that already very small external per- 
turbations can have a significant impact on effective clus- 
ter states. Typically, the effective multi-site interactions 
yielding the effective cluster state arise in a high order in 
perturbation theory (here order 4). As a consequence, 
any external perturbation acting in the effective low- 
energy model in a lower order (here order 1) represents 
a strong constraint for the effective implementation of a 
1-WQC. This effect is present on any lattice and in any 
dimension for the problem studied in this work. 

The physical mechanism leading to the dramatic loss 
of fidelity is actually very different for the two pertur- 
bations we have considered. The external Z-field leads 
to a polarization of the ground state and therefore to 
a reduction of entanglement. The additional Ising cou- 
pling induce thermal and quantum fluctuations due to a 
quantum phase transition. 

The qualitative aspects of our work are relevant for a 
much broader class of problems: any effective low-energy 
model which is derived perturbatively and which con- 
tains dominant multi-site interactions is expected to be 
affected by external perturbations. The physical reason is 



that effective n-site interactions arise typically in order n 
while it is likely that external perturbations exist which 
act non-trivially on the effective low-energy degrees of 
freedom already in a lower order. 

A prominent example is Kitaev's honeycomb model 
which contains the so-called toric code as an effective 
low-energy model perturbatively in order 4 [15]. The 
toric code is a topological stabilizer code consisting solely 
of 4-spin interactions. One can easily show that exactly 
the same kind of external perturbations studied in this 
work give again rise to operators in the effective model 
in order one perturbation theory causing a breakdown of 
the topological phase for small external perturbations. 

In the light of the severe constraints found in this work 
for the realization of effective cluster states, let us finally 
mention concepts for MBQC using elementary entities 
with larger spins which represents a promising route for 
future research [16-19]. 
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SUPPLEMENTAL MATERIAL 



Spectral properties of H 



The spectrum of H (and H diag ) is the product of the spectra of the local H l ° c terms. In Tab. I one can see the 
eigenenergies of the H l ° c eigenvectors |0) , |15) . The unique ground state 



IV-. 



|0>„ 



of H dmg is characterized by all lattice sites being in the |0) state. The first excited space is iV-fold degenerated and 
is spanned by the states (/i G C) 



(a0\ 

EX/ 



>„ 1°)- 

uec/{fi} 



For the energy gap AE = E\^ 
one therefore finds 



Ejo) between the ground state and the (A 7 "- fold degenerated) first excited space 




AE 



It is independent of the system size. From the Taylor expansion of AE one sees, that it is of order 4 in A a 



AE: 



5 / X X z 



32 V £ 



21 
512 




All of the first TV excited spaces can be characterized by some lattice sites being in the |1) state and all other 

lattice sites being in the |0) state. They form the space S of the effective qubits. The n-th excited space is (^) times 
degenerated and the gap between the n-th and the (n + l)-th space is AE. All high-energy states are characterized 
by at least one lattice site being in one of the states |2) , |15). In Fig. 4 you can see the lower part of the spectrum 
for the example N = 4. 



State 


E \i) u 




-2^2 + 2^+2^ + 1 


ID, 


-2g-2g ^Jl + ^r 


\ 2 l 


-2g ^2 + 2^-2^ + 1 


|3>„ , |4>„ 


-2g Kz 


|5>„ 


2g-2g^Jl + ^ 


m u ,\7) u ,\8) u ,\9) u 





|io> M 


-2g + 2g <Jl + ^ 


|11>„ , |12>„ 


2g \ xz 


|13> M 


2^2 + 2^-2^ + 1 


|14> M 


2g + 2g^l + ^ 


|15> M 


2^2 + 2^ + 2^ + 1 



Table I. Eigenenergies of the H l ° c eigenvectors |0)^ , |15) 



6 




Figure 4. Lower part of the spectrum of H for N = 4 



Thermal fidelity of H z 



To calculate the fidelity per site for finite temperatures for the case of an external Z-field one must note, that 
the points (b) and (c) are no longer valid to justify the approximation, since the ground-state fidelity quickly de- 
creases with h z and since K is no longer a conserved quantity Because of (a) it still is possible to approximate 
e -P E \*z) \ (i z \R z \+) | 2 w Vz > 2 and Z((3) « e~^ E \ Q z) + e - ^ 1 *), where is the 16 x 16 Matrix of the 
eigenvectors |0 2 ) , 1 15^) . One finds 

d - 1+e -^ g j <wi+> i 2 + rri^i <i.i^i+> i 2 - w 

The energy gap AE Z = E\ lz ^ — -E'|o z >? the fidelity d = \ (1 Z \R Z \+) | 2 of the first excited state , and the fidelity per 
site d for the temperatures T = 0.0001 g/ks and T = O.OOlg/ks are shown in Fig. 5. By comparing with the fidelity 
per site for T = (shown in the main paper) one sees, that for finite h z thermal fluctuations play a minor role, as the 
gap is increased by the field. 



Perturbation method by Takahashi 

To calculate the ground- state fidelity per site of the TFIM on a square lattice with the polarized state, the fidelity 
per site of H zz with the logical cluster state, and the energy gap of H zz as series expansions, we used a perturbation 
method introduced by Takahashi [1]. In terms of this method a Hamiltonian H = Hq + V is transferred by a base 
transformation T into an effective Hamiltonian i^ e ff, so that it only acts in one of the Hq eigenspaces, while the 
spectral properties are preserved. 

Let P be the projector into a H eigenspace L with eigenvalue and P be the projector into the according 
perturbed H eigenspace L. The Operator T is then defined by 

T := PP(PPP)-^ 

with r^T = P. The projector P can be expressed as 

P = P _Y^ J2 S kl VS k2 V...VS k "+ 1 (10) 

n=l k\+k2+----\-k n j t i=n, ki>0 
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Figure 5. Energy gap AE Z — E\ lz ) — E\ 0z ) (top left), fidelity per site d — \ (1 Z \R Z \+) | 2 of the first excited state (top right), 
and fidelity per site d (Eq. 9) for the temperatures T = 0.0001 g/ks and T — 0.001 g/ks (bottom left and right) for the case 
H z of a perturbation with an external Z-field. 
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with S° := -P and S k := (^wrjf^J ■ Together with 

(pppyi = p+jr (2 ? 2 ~J, )! W - p)p) n (ii) 

n—l ^ '" 

it is now possible to express Y order by order as a power series of V. Up to order 3 the coefficients of Y = J^rW are 
given by 

r co) = p 5 

r (1) = svp, 
r (2) = svsvp - s 2 vpvp - -pvs 2 vp, 

2 

r (3) = s 3 VPVPVP - S 2 VSVPVP - SVS 2 VPVP 
+ ^PVS 3 VPVP - S 2 VPVSVP + SV SVSVP 

--PVS 2 VSVP - -SVPVS 2 VP 
2 2 

--PVSVS 2 VP + ^PVPVS 3 VP. 

Within the context of this work we calculated the coefficients up to order 13 (7030571 coefficients). 
We now transform the Schrodinger equation into the effective basis 

H\^ i )=E i \^ i ) 

^ eff r t |^> = EiT* |^> 

with H e fi := Y^HY. Every eigenvector of H in L therefore corresponds to an eigenvector of H e R in L with the 
same eigenvalue. It is possible to calculate the perturbed eigenvector |^) with the base transformation Y from the 
according H e ^ eigenvector Y^ |^). 

To express the effective Hamiltonian H e ^ = as a power series of V we exploit that [H, P] = 

#eff = (PPP)-^PPHPP(PPP)-^ 

= (PPP)-ip(Ho + V)PP(PPP)-^ 
= P + (PPP)-* PVPP(PPP)-^ . 
Together with Eqs. 10-11 one finds up to order 4 

rr(0) _ -(0) p 

H$ = PVSVP, 
= PV SVSVP - ^PVPVS 2 VP - ^PVS 2 VPVP, 

= PVSV SVSVP + ^PVPVPVS 3 VP - ^PVPVSVS 2 VP 

--PVPVS 2 VSVP - -PVSVPVS 2 VP - -PVSVS 2 VPVP 
2 2 2 

-^PVS 2 VPVSVP - ^PVS 2 VSVPVP + ^PVS 3 VPVPVP. 

Within the context of this work we calculated the coefficients up to order 14 (5394321 coefficients) and for the 
special case PVP = up to order 17 (2490673 coefficients). 
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Application of Takahashi's method to the energy gap 

We applied this perturbation method to calculate the energy gap between the ground state of H zz and its first 
excited state up to order 5. We start from the representation of H in the basis H diag and we define 

H := # diag , 

A method to calculate series expansions for properties of excited states is based on a work by Gelfand [2]. The 
idea is to interpret excitations on the lattice as quasi-particles and to describe their dynamics. In this quasi-particle 
picture, the ground state \^h) of H dmg corresponds to the vacuum state and the state |?/4x) corresponds to a state 
with a single particle at lattice site /i = (x, y) where (x, y) are Cartesian coordinates. By switching on the perturbation 
part, dynamics between quasi-particles are induced. Using the effective form H e ^ of H zz acting solely in the space of 
the first excited states, we define the hopping elements 

== {^ xMv) \H eS \^ b) } a,b,x,y€ Z. 

To this end, it is important that the perturbation V is a sum of local operators each linking two neareast-neighbor 
lattice sites. A linked process is then defined as a linked operator sequence of these local operators. The linked-cluster 
theorem [3, 4] states that only linked processes contribute to the one-particle amplitudes. It is therefore possible 
to retrieve the hopping elements in a given order (in our case order 5) in the thermodynamic limit, although the 
calculations are limited to finite systems. 

To achieve this the finite systems have to be constructed large enough to contain all linked processes of up to 5 
bonds. It also follows, that \x\ + \y\ > 5 => t XjV = 0, as there is no order 5 process linking the lattice sites (a, b) and 
(a + b + y) for this case. 

The calculation can further be simplified: 

1. For symmetry reasons it is 

^x,y t—x,y t x —y t— x —y , 

so only one calculation is necessary for such hopping elements. 

2. An operator Z rn ^Z rn ^ 2 ) of V anti-commutes with K m m and K m (2). It follows that the application of such an 
operator inverts both eigenvalues of the two neighboring lattice sites. To contribute to a hopping element, 
each lattice site must therefore be touched an even number of times by such operators. The only exception to 
this rule are the lattice sites (a, b) and (a + x, b + y) for (x,y) ^ (0,0). They must be touched by an uneven 
number of operators. This allows to further decrease the size of the finite systems. Examples for some hopping 
elements can be found in Fig. 6. 

3. If h is an operator of the power series representation of H e R one can calculate its contribution to a hopping 
element by calculating h |^ex^) an d multiplying the result with (^ex^'^^I- But to save ^ me an d memory it 
is useful to split the coefficient h = h\h r into two parts with roughly the same number of V Operators. The 
contribution is then calculated by multiplying the two vectors h r |^ex^) an d (V4x" x,b+y " > | hi. 

From the hopping elements we calculate the dispersion of the first excited mode with a Fourier transformation. For 
the square lattice one finds 

u(k x ,k y ) =t 0l o- E + ^2 t XjV cos(k x x + k y y), (12) 

where E Q is the ground-state energy of the finite lattice used to calculate to,o- The dispersions for the case X xz = 0.5 g 
and for different X zz are shown in Fig. 7. The minimal excitation energy (the gap) is at k x = k y = tt for negative X zz 
and at k x = k y = for positive X zz . The coefficients of the series expansion of the gap up to order 5 are shown in 
Tab. II for different X xz . The first-order coefficients equal the coefficient for the TFIM gap. The absolute value of all 
other coefficients are decreased for finite X xz . The actual value of the critical A^lcrit is therefore slightly larger then 
the critical point calculated using the low-energy approximation (see Fig. 8). 
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Figure 6. Examples for finite systems to calculate one-particle hopping elements t x , y from the lattice site (a, b) to (a + x, b+y) or 
vice versa are shown. These two lattice sites are highlighted by a green (gray) background. The finite systems are constructed 
such that all connected graphs of up to five bonds (a single bond can be counted multiple times) are contained, where each 
lattice site is touched by an even number of bonds. The only exception to this rule are the lattice sites (a, b) and (a + x, b + y) 
which must be touched by an odd number of bonds for (x, y) 7^ (0, 0). The hopping elements of this systems equal the hopping 
elements of an infinite system up to order 5. 
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Table II. Coefficients of the energy gap AE + a ^f- + c 2 (^) 2 + c 3 (^) 3 + c 4 (^) 4 + c 5 (^) 5 of H zz up to order 5. The first 
row shows the coefficients for the TFIM (limit X xz 0). It is fi := {y^zY ■ 



Application of Takahashi's method to the ground-state fidelity 



We have applied the perturbation method to calculate the ground-state fidelity per site of H zz with the logical 
cluster state up to order 4 and the ground-state fidelity per site of a TFIM on a square lattice with the polarized state 
up to order 12. For H zz , we again use the basis £ diag and define 

H := # diag , 
V := —X zz R^cZm(l)Zm(2)RC' 

For the TFIM, we use 

fi£C 



y '■= ~^zzC 2^ ^m(l)^m(2)- 
meM 

The ground state \^h) = °^ * s m both cases unique and therefore also the only eigenvector of the 

according H e R. The wave functions of the perturbed ground states can therefore be calculated by applying the base 
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Figure 8. Critical \ zz \cvit in dependence of X xz calculated with the low-energy approximation (black solid line) and corrected 
by high-energy contributions (red dashed line). The high-energy corrections are calculated by a dlogPade [2,2] approximation 
of the order 5 series of the energy gap. 

transformation T to |0) M - For the ground-state fidelity per site one finds 

d Hzz (\^cs) , WhJ) = Jim Vl (^cs aS |r|0) £ | 2 , 

d TF iM(|0) £ , |^tfim» = Jim^ )/| (0| £ r|0) £ | 2 . 

Again it is possible to exploit the linked-cluster theorem to limit the calculation to finite systems. For example it is 
sufficient to calculate on a 5 x 5 lattice with periodic boundary conditions (system lies on a torus) for the order 4 
calculation, or a 13 x 13 lattice for the order 12 calculation respectively. 
The calculation can further be simplified. 

1. Analog to the second point of the last section only processes can contribute, where each lattice site is touched 
by an even number of operators located on the bonds. It is easy to see, that this is only possible if the overall 
number of bonds is even. For the calculations only coefficients of even order must be considered. 

2. The 5x5 lattice necessary for the order 4 calculation can be reduced to a 3 x 5 lattice by bounding it periodically 
into a "brick wall" structure (see Fig. 9). 

3. To reduce the size of the 13 x 13 lattice for the order 12 calculation another technique was used. The coupling 
parameter X zz was replaced by the parameters A| and A_, so that vertical bonds are associated with A| and 
horizontal bonds are associated with A_. Using the rotational symmetry of the lattice the calculation of the 
order 12 coefficient c\2 is reduced to the calculation of the A| and A_ coefficients 

c 12 XH = c& 12) (Af 2 + X 1 - 2 ) + c£ 10) (Af \]_° + A[°A 2 _) + <& 8) (Af A 8 _ + Af A 4 _) + cg- 6 >Af A 6 _. 

For the separate calculations of the coefficients c^' 12 \ c^' 10 ^ c^' 8 ^ and c^' 6 ^ periodically bounded systems of 
size lx 13, 3x11, 5x9 and 7x7 are sufficient. 
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Figure 9. Periodic lattice with "brick wall" boundary conditions. All linked processes of at most 4 bonds cannot form closed 
loops crossing two opposed lattice boundaries. The order 4 fidelity per site of this lattice is therefore identically to the fidelity 
in the thermodynamic limit. 



4. It is (0\ c S = 0. For the calculation of (Ol^riO)^ all coefficients with an S'-operator at the left side can therefore 
be left out. 

5. Analog to point 3 of the last section it is possible to reduce the calculation of the contribution (0|/:7|0)£ of 
a coefficient 7 = 7^7 r of the power series of T to the calculation of (0^7^ and j r \0} c and their product. In 
principle this can be also applied to the calculation of (i/>Qs S \'y\0) £• Due to the size of the representation of 
(ipcs g \ m tne basis H dia s this does not save time or memory. It is better to calculate 7 (0)^ and multiply it with 
<^ s as | directly. 

6. For symmetry reasons the contribution of a process only depends on the shape of the process and not on its 
position on the lattice. It is therefore possible to limit one of the ^/-operators (advisably the first) to a single 
and fixed bond. The resulting value is then the fidelity of the system divided by the number of bonds. 

The resulting ground-state fidelity per site of the TFIM on a square lattice with the polarized state is given by 



^TFIM 
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812949139 
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The coefficients of the ground-state fidelity per site of H zz with the logical cluster state are shown in Tab. III. One 
sees, that the absolute values of the coefficients are decreased for finite X xz . The actual ground-state fidelity per site 
of H zz is therefore slightly larger then the fidelity calculated with the low-energy approximation. The fidelity g?tfim 
and the fidelity of H zz for different \ xz are plotted in Fig. 10. 
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Table III. Coefficients of the ground-state fidelity per site d Hzz = \ (0\R\+) | 2 + c 2 (^) 2 + c 4 (^) 4 of H zz with the logical 
cluster state for different X xz . The first row shows the coefficients for the TFIM (limit X xz — > 0). 




Figure 10. Ground-state fidelity per site d of H zz calculated with the low-energy approximation (cZtfim, shown as solid line) 
and corrected by high-energy contributions for different X xz (dashed and dotted lines). The low-energy fidelity is calculated up 
to order 12 and the corrections are calculated up to order 4. 



